12/13/05 TUE 11:33 FAX 5152881338 



MCKEE VOORHEES & SEASE 



PTO 



Application of Mathematical Tools for 
Metabolic Design of Microbial 
Ethanol Production 



@022 



Vassily Hatzirnanikatis,* Marcel Emmerling, Uwe Sailer, James E. Bailey 

Institute of Biotechnology. ETH Zurich, CH-8093 Zurich, Switzerland; 
telephone: 411 633 31 70; fax: 411 633 10 51 

Received 7 July 1997; accepted 30 Jufy 1997 



Abstract: Many attempts to engineer cellular metabolism 
have failed due to the complexity of cellular functions. 
Mathematical and computational methods are needed 
that can organize the available experimental information, 
and provide insight and guidance for successful meta- 
bolic engineering. Two such methods are reviewed here. 
Both methods employ a (log)linear kinetic model of me- 
tabolism that is constructed baaed on enzyme kinetics 
characteristics. The first method allows the description of 
the dynamic responses of metabolic systems subject to 
spatiotempore! variations in their parameters. The sec- 
ond method considers the product-oriented, constrained 
optimization of metabolic reaction networks using 
mixed-integer linear programming methods. The optimi- 
zation framework is used in order to identify the combi- 
nations of the metabolic characteristics of the glycolytic 
enzymes from yeast and bacteria that will maximize 
ethanol production. The methods are also applied to the 
design of microbial ethanol production metabolism* The 
results of the calculations are in qualitative agreement 
with experimental data presented here. Experiments and 
calculations suggest that in resting Escherichia coli cells, 
ethanol production and glucose uptake rates can be in- 
creased by 30% and 20%, respectively, by overexpres- 
slon of a deregulated pyruvate kinase, while increase in 
phosphofmctokinase expression levels has no effect on 
ethanol production and glucose uptake rates. ©1998 John 
Wiley & Sons. Inc. Biotechno! Bioeng 58: 154-161. 1998. 
Keywords: central carbon pathways; metabolic optimiza- 
tion; ethanol production 



INTRODUCTION 

Current knowledge of biochemical systems is composed of 
a vast set of data that accumulates at an increasing rate. 
Advances in analytical methods and development of sophis- 
ticated instrumentation have provided the tools that allow us 
to know more than we can understand. Comprehension is 
complicated by the high complexity of living organisms. 
This complexity increases from unicellular organisms to 
multicellular structures, such as tissues and organs, to dif- 
ferentiated organisms. 
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Recombinant DNA and other modern; genetic technology 
allow the introduction of precise changes in the genetic 
determinants of cell function and, thus, the engineering of 
metabolic activities for novel and/or improved functions. 
The application of molecular genetics and engineering tools 
for the useful manipulation of cellular processes is defined 
as metabolic engineering (Bailey, 1991). This definition 
identifies the engineering virtues of the manipulation of 
metabolism: putting together the available information, em- 
ploying intuition, and tackling the problem. However, as in 
many engineering approaches, solution of the problem has 
not been always achieved. Moreover, the discovery of some 
of the most profitable bio technological applications has 
been serendipitous- 
Current approaches and methodologies in the biological 
sciences, while providing understanding of isolated cellular 
processes and subsystems, do not promote understanding of 
the simultaneous contributions of these subsystems to over- 
all cellular metabolism. It is the organization of these sub- 
systems, which are themselves relatively small but elaborate 
networks, in large complex networks which cause intuition 
and piece-wise knowledge to fail as tools for rational meta- 
bolic engineering (and, similarly, for accomplishing mo- 
lecular medicine). Moreover, the high dimensionality of 
these networks offers an almost infinite number of alterna- 
tive approaches toward the achievement of a goal. Most of 
the metabolic networks considered for engineered manipu- 
lation are composed of enzyme-catalyzed reactions with 
complicated stoichiornetry, nonlinear kinetics, and superim- 
posed regulatory structures. These regulatory structures in- 
clude interactions of enzymes with compounds in the reac- 
tion networks, which often are not reacxants or products of 
the reaction catalyzed by the enzymes in question, but serve 
to modulate the catalytic activity of the enzymes. Manipu- 
lation of a metabolic subsystem will change the metabolite 
concentrations and the metabolic activities of the sub- 
system- This manipulation and the resultant subsystem 
changes often drive responses in the subsystem's cellular 
environment that have not been taken into consideration. 

The development of tools and frameworks that will or- 
ganize the available biological knowledge and that will help 
in the analysis and design of metabolic networks is of im- 
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mediate importance. These tools should be able to screen 
efficiently through an almost unlimited set of cellular modi- 
fications (realizable by genetic engineering technology), 
and report a small set of most promising options that can be 
further tested for their potential in the development of new 
products and improved bioprocesses. 

MATHEMATICAL METHODS FOR ANALYSIS AND 
DESIGN OF METABOLIC SYSTEMS 

Mathematical models of cellular processes have been used 
successfully in the past for the analysis of these processes, 
and they have provided useful initial directions for genetic 
improvements of the process of interest The realization that 
mathematical models are the only way that net conse- 
. quences of simultaneous, coupled, and often counteracting 
processes can be evaluated consistently and quantitatively 
has led to the growth of mathematical modeling in many 
biological and biotechnological areas. The analysis of these 
models requires development of associated mathematical 
and computational methods. Two mathematical and com- 
putational methods are proposed here as tools for obtaining 
useful guidance for experimental application of metabolic 
engineering. 

One of the most common targets of metabolic engineer- 
ing is the manipulation of metabolic functions, such as 
metabolic reaction rates and intracellular metabolite con- 
centrations, in desired directions. Prediction of the changes 
in metabolic functions, which will occur after a change in 
any parameter affecting the metabolism, requires a kinetic 
representation relating fluxes to metabolite concentrations 
and metabolic parameters. In general , such models are not 
available, and the kinetic model must, therefore, be devel- 
oped from measurements of reaction rates and metabolite 
concentrations. A variety of formulations are possible. A 
(log)linear kinetic model for metabolic reaction networks is 
presented here (Hatzimanikatis and Bailey, 1997a). The 
model does not require detailed information about the ki- 
netic mechanisms of the reactions. It simply employs 
knowledge about the strength of interaction among the rates 
of enzyme-catalyzed reactions and the concentrations of the 
various metabolites, substrates, and regulators of the meta- 
bolic network. This kind of information can be derived from 
experimental and theoretical methods developed around 
metabolic control analysis (MCA), a sensitivity analysis 
framework that allows determination of quantitative indices 
such as control coefficients and elasticities (Cornish- 
Bowden and Cardenas, 1990; Fell, 1992; Kacser and Bums, 
1973; Schlosser et aL 1993). Comparative studies between 
(log)linear models and nonlinear models based on common 
reaction mechanisms have been performed, demonstrating 
the satisfactory accuracy of (log)linear models in approxi- 
mating the dynamic responses of metabolic networks to 
changes in metabolic parameters (Hatzimanikatis and Bai- 
ley, 1997a). 

Unsuccessful attempts to engineer cellular metabolism by 
simply manipulating the amount of various enzymes sug- 



HATZIMANIKATIS ET AL: MATHEMATICAL TOOLS FOH 



gest that engineering of the regulatory characteristics of the 
enzymes in a metabolic network offers greater potential for 
the achievement of desired metabolic properties. In many 
examples of small reaction networks wilfa simple regulatory 
structures, changes in the regulation of one or two enzymes 
improved product formation significantly (Bailey, 1991; 
Katsumata and Ikeda, 1993). m large metabolic networks, 
or those closer to central carbon metabolism, the regulatory 
structures tend to be more elaborate. Moreover* the experi- 
mental difficulties in modifying regulatory interactions add 
another degree of complexity, because: the common trial- 
and-error experimental approach is unfeasible. A produc- 
tion-oriented optimization framework has been recently de- 
veloped (Hatzimanikatis et aL, 1996). the framework con- 
siders the optimization of the performance of a metabolic 
pathway with respect to changes in the amounts of enzymes 
in the pathway and in modifications in the regulatory char- 
acteristics of those enzymes. The optimization study is un- 
dertaken using the kinetic description provided by the (lo- 
g)linear kinetic model, used in this case to define constraints 
on the optimization in the form of steady-state mass bal- 
ances of intracellular metabolites. Consideration of changes 
in enzyme regulation which abolish initially present inhibi- 
tion or activation introduces binary decision variables into 
the optimization, resulting in a mixed-integer programming 
problem. 

The development of these mathematical methods has 
been extensively presented elsewhere (Hatzimanikatis et al. a 
1996; Hatzimanikatis and Bailey, 1997a). We will present 
here an application Co microbial ethanol iproduction, and the 
results of this analysis will be compared to experimental 
data. 



METABOLIC ENGINEERING OF MICROBIAL 
ETHANOL PRODUCTION 

Ethanol production from sugar mixtures; by microorganisms 
has been the subject of extensive research (Ingram and 
Doran, 1995). Saccfiaromyces cerevlsiae and Zymomonas 
mobilis have been the most widely used. organisms for etha- 
nol production. However, in recent years Esclierichia coli 
has been genetically engineered for ethanol production 
(Ohta et a!., 1991). The main pathway involved in ethanol 
production is the glycolytic pathway; Figure 1 Shows this 
pathway and its regulation in S. cerevisiae. The complexity 
of glycolysis and its regulation has been the main obstacle 
in many experimental attempts to increase ethanol produc- 
tion and, in general, to manipulate its functions by meta- 
bolic engineering (Boles et aL, 1993; Sbhaaff et aL, 1989). 
Mathematical modeling of the glycolytic pathway in various 
organisms and computer simulations have also been used 
for a better understanding of the performance of the path- 
way and its responses to genetic manipulations (Gallazo and 
Bailey, 1990; Hatzimanikatis and Bailey, 1997b; Heinrich 
et aL, 1977; Joshi and Palsson, 1989; Schlosser et aL, 1994). 

JC DESIGN OF MICROBIAL ETHANOL PftOfiuCTION 1 55 
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Figure i. Anaerobic fermen ration pathway of nongrowiag yeast> 5flC- 
chttromycex cerevisieu. with glucose as the sole carbon source. Chemical 
species: intracellular glucose; G6P, glucose 6-phosphate; l*6P, fructose 
6-phosphatc: FdP, fructose 1 ^-diphosphate: 3PG, 3-phosphogly cerate; 
PEP, phosphoenolpyruvate; POL, polysaccharides; GRO, glycerol: £TOH, 
ethanol: ATP, adenosine triphosphate; ADP. adenosine diphosphate; AMP. 
adenosine monophosphate. Pathway steps and enzymes; In, glucose up- 
take; HK, hexokinase; equilibrium step; PFK, phosphofructokinase; 
GAPD, glycerol dchydc 3-phosphate dehydrogenase; K u equilibrium step; 
PYK, pyruvate kinase; GRO, glycerol production; POL» polysaccharide 
production: ATPase, net ATP consumption: AK, adenylate kinase. Solid 
arrows indicate reaction steps, dashed arrows indicate activation and dotted 
arrows indicate inhibition. 



A (Log)linear Modal of Microbial 
Ethanol Production 

A mathematical model for the description of ethanol pro- 
duction by nongrowing 5- cerevisiac has been constructed, 
and found to be consistent with experimental observations 
(Gallazo and Bailey, 1990; Schlosser et aU 1994). This 
model uses nonlinear rate expressions based on kinetic pa- 
rameters estimated from in vitro studies, as well as expres- 
sions parametrized by in vivo perturbation experiments. The 
pathway features seven regulatory loops of activation and 
inhibition, and the kinetic expressions are strongly nonlin- 
ear. Based on this nonlinear model, we constructed a (lo- 
g)lincar model around a reference steady state. The corre- 
sponding parameters of the (log)lmear model, i.e., the stoi- 
chiometric and elasticities matrices* are presented in the 
Appendix. 

We performed a comparison between the (log)lmear and 



the nonlinear models with respect to; their dynamic re- 
sponses to changes in metabolic parameters. Characteristic 
results for the responses of the two models to step and 
pulscd-periodic variations in glucose uptake are presented 
in Figure 2, and they demonstrate the excellent accuracy of 
the (log)lincar model The final steady-state differences for 
step responses are due to inherent limitations in MCA cal- 
culations of new steady states. , 

Excellent agreement between response characteristics of 
the (log)linear model and of the corresponding nonlinear 
model has also been observed for variotis complicated path- 
ways (Hatzimanikatis and Bailey. 1991a). This agreement 
suggests that, in the absence of a nonlinqar model, and when 
MCA quantities can be experimentally idetermined or esti- 
mated, die (log)linear model can be used as a good first 
approximation for the analysis of the- dynamic response 
characteristics of metabolic systems. 

Metabolic systems, like any nonlinear dynamic system, 
can exhibit a shift in time-averaged properties (such as 
metabolic fluxes and control coefficients) in response to 
fluctuations with zero-time-average value. Such fl actua- 
tions can arise, for example, in bioreactors when mixing is 
nonideal, a situation common in large-scale systems. Based 
on the excellent approximation provided by the (log)linear 
model, a method has been developed for the estimation of 
the performance of metabolic systems i subject to periodic 
spatio temporal variations of the system parameters and pro- 
cess operating conditions. Application of the method to the 
ethanol production pathway has shown that, for periodic 
fluctuations in extracellular glucose concentration, time- 
average flux control coefficients vary strongly, and not 
monotonically, with the period of the external fluctuations 
(Hatzimanikatis and Bailey, 1997a). These observations are 
particularly significant because they imply that decisions 
taken at a certain stage of process scale-up, with respect to 
which enzyme in a pathway should be genetically manipu- 
lated, could have no effect, or even be counterproductive at 
another scale. 
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figure 2. Dynamic responses of the ethanol specific production rate of the yeast glycolytic pathway. Solid lines correspond to* the nonlinear model and 
dashed lines correspond to the (logjlincar model. A & C. Responses to 5t<?p changes of the glucose uptake. J; +20%; II: + J09b; Tit -10%; IV: -20% (from 
Hataraanikalis and Bailey, 19971>.) B. Responses to periodic variation of the input flux with period T = 1 and v^CO = 1.2v 4fc « for nT « t < (n + 0.5)T 
and v*, CO = 0.8 ^ # for [n + 0.5)r *r < {n + 1)7'. j 
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Optimization of Microbial Ethanol Production 

The kinetic description provided by the (log)linear meta- 
bolic model has been used for identification of the changes 
in the amounts of the emymes of the glycolytic pathway 
and of the modifications in the regulatory characteristics of 
those enzymes that will maximize ethanol production. 

The problem of maximizing microbial ethanol production 
has been addressed in two different contexts: the analysis 
problem and the synthesis problem. In the context of the 
regulatory structure present in yeast, which can be modified 
by inactivating any of the regulatory loops, the analysis 
problem is postulated as: 

Which of (he existent regulatory loops should be inactivated, and 
what associated changes should be made in enzyme expression 
levels, in order to maximize Specific ethanol production? 

In a more extended context, a regulatory superstructure is 
postulated which combines the regulatory characteristics of 
the glycolytic pathway in the bacterium, E. coli and yeast, 
and the synthesis problem considers the following ques- 
tions: 

flow should the regulatory loops present In both E. coli and yeast 
be combined, and what associated changes should be tmtdc in 
enzyme expression levels in order to maximize specific ethanol 
production? 

The number of the regulatory loops in the synthesis prob- 
lem is equal to 7 and tbe number of the enzymes that can be 
manipulated is equal to 8. Therefore, there are 2 (74 * } = 
32,768 alternative solution structures. In the synthesis prob- 
lem, the number of the regulatory loops is equal to 10 and 
the number of the enzymes that can be manipulated is equal 
to 8, resulting in 2 fl04 " 8) — 262,144 alternative solution 
structures. The highly combinatorial nature of the problem 
manifests the need for a mathematical and computational 
method that can screen for the most promising initial strat- 
egics. 

The optimization calculation for the analysis problem in- 
dicated that only two activation loops should be retained, 
and three enzymes should be overexpressed in order to 
maximize ethanol production, giving an increase of 100% 
relative to the reference State. The optimal regulatory struc- 
ture and enzyme manipulation pattern is illustrated in Fig- 
ure 3. 

Solution of the synthesis problem, which begins from a 
pathway containing the superimposed enzymatic regulation 
of E. coli and of S. cerevisiae (Fig. 4) suggested the regu- 
latory structure presented in Figure 5 and overexpression of 
three enzymes, with a higher overexpression level for py- 
ruvate kinase (Pyk), that could increase ethanol production 
by 114%. In the optimal glycolytic pathway suggested by 
the calculations, phosphofructokinase (Pfk) features regula- 
tory characteristics from £. coli, and Pyk is affected by only 
one regulatory loop, activation by FdP, common in both E. 
coli and yeast. 

We could, therefore, draw the conclusion that, in an etha- 
nol-producing E. coli strain, overexpression of a Pyk that 
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Figure 3. The optimal regulatory structure and enzyme manipulation 
calculated for the analysis problem for optimization of yeast ethanol pro- 
duction. Thiclc solid arrows Indicate increases in enzyme levels for the 
respective reaction. Abbreviations for chemical species and pathway steps 
are the some as in Figure 1. 

lacks inhibition from ATP might increase the ethanol pro- 
duction rate. Moreover, the optimization solution suggested 
that overexpression of Pfk will not increase ethanol produc- 
tion significantly. These conclusions are strictly qualitative 
because the main parameters used. i.e.. substrate elasticities, 
are based on a kinetic model for yeast, cand certain stoichi- 
ometric structures from yeast, such as glucose phosphory- 
lation and glycerol production, which are not present in E. 
coli t 

The suggested solution for die synthesis problem was 
compared with results from laboratory experiments with an 
ethanol-producing E. coli strain under anaerobic, nongrow- 
ing conditions. The ethanol-producing JS. coli strain KO20 
(Ohta et al. t 1991). obtained from Prof.; L.O. Ingram (Uni- 
versity of Florida), was used as the host in these experi- 
ments. The strain is characterized by ond chromosomal copy 
of the artificial pet operon encoding the Zymomonas mobilis 
genes for alcohol dehydrogenase II andjpyruvate decarbox- 
ylase, integrated within the pyruvate formate lyase gene. 
Thus, under anaerobic conditions, this strain produces 
mainly ethanol as the catabobc product of glucose utiliza- 
tion. Based on the £ coli expression vector pTrc99A. (Phar- 
macia LKB, Uppsala, Sweden), several plasmids carrying 
homologous and heterologous pfk and pyk genes under the 
inducible tac promoter were constructed, allowing moderate 
overexpression of the corresponding enzymes (2-10 times). 
The enzymes expressed in combinations were the E. coli 
PfkA and PykF, and the heterologous iion-ailosteric, non- 
regulated Pfk from Lactobacillus bulgaricus (Branny et al., 
1993; Le Bras et al., 1991) and Pyk from Bacillus stearo- 
thermophilus (Sakai and Ohta, 1993), which lacks any ac- 
tivation by FdP, but is allosterically activated by AMP and 
ribose-5-phosphate (Sakai et al., 1986). Details of the ex- 
perimental procedure can be found elsewhere (Emmerling 
et al., in preparation). In all experiments reported here, cells 
were cultivated under anaerobic, non-growmg conditions. 
In order to determine specific glucose .uptake and ethanol 
production rates, the extracellular glucose and ethanol con- 
centrations were determined at successive time points. 

Specific ethanol production rate was increased by 1 1% in 
the strain overexpressing the homologous pyk (PYKec), and 
by 17% in the strain overexpressing tike heterologous pyk 
from B. stearothermophilus (PYKbs), under conditions of 
low initial glucose concentration and: for cells harvested 
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Figure 4. The glycolytic pathway with the regulatory superstructure composed by superimposing regulatory characteristics from yeast (dashed arrows) 
and from E- coli (dotted arrows). Dashed^dotied arrows indicated regulatory loops common to the two orgasms. Plus signs indicate activation and minus 
signs indicate inhibition for the respective loops. Abbreviations for chemical species and pathway steps are! the same as in Fig' 
ure 1. 



from mid-exponential growth phase (Fig. 6A). Moreover, 
specific glucose uptake rate was increased by 10% in strain 
FYKbs compared to the identical uptake rates of the control 
strain carrying only the pJasmid and the PYKec strain (Fig. 
6B). UndeT these conditions, constant glucose uptake rate 
and ethanol production rate were observed over the time 
course of the experiment (1 h) for all strains, while a de- 
crease in the glucose uptake rate over lime was observed for 
all strains under conditions of high initial glucose concen- 
tration, and for cells harvested from late exponential growth 
phase (Figs. 7A & B). However, in the PYKbs strain, glu- 
cose uptake rate was higher by 18% as compared to the 
control and the PYKec strains (Fig. 7B). Tn addition, no 
apparent decrease in the specific ethanol production rate 



was observed in the PYKbs strain, contrary to the control 
and the PYKec strains in which ethanol production rate 
decreased over time (Fig. 7A), Thus, under the conditions 
investigated, moderate overexpression :of the heterologous 
pyk from B. stearothermophilus increases the glycolytic 
flux in E. coli engineered for ethanol production by 10% to 
20%» and the final ethanol yield by 20& to 30%. 

These presented experimental results are in qualitative 
agreement with the suggested solution! from the synthesis 
problem. Strain PYKbs has a regulatory structure similar to 
the calculated optimal structure. However, Pyk from B. 
stearotfiermophilus is not activated by FdP, whereas the 
optimal structure suggests the desirability of this interac- 
tion. 
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Figure 5. The optimal regulatory structure and enzyme manipulation calculated for the synthesis problem for maximizing ethbiOl production. Dashed- 
dotted arrows indicate regulatory loops common to the two organisms. Plus signs indicate activation and minus signs indicate inhibition for the respective 
loops. Abbreviations for chemical specie? and pathway steps arc the same as in Figure 1. 
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Ficurc 6. Specific glucose uptake (A) and ethanol production rate (B) at 
0-5% initial glucose concentration of anaerobic, resting £ coti KO20 
strains with the plasmidp7>c9 a A only (control), ove^cxp^essingp/l/ , from 
£ coti (PYKec). pnd overexposing pyk from B. sUarothennophiius 
(PYKbs). Cells were harvested from anaerobic mid-exponcntiaJ growth 
phajc cultures. The indicated values represent the mean value of at least 
two independent experiments, and error bar? indicate the standard devia- 
tion. 



CONCLUSIONS 

A precise quantitative description of metabolism will, in 
most cases, be impossible, at least in the foreseeable future. 
However, careful mathematical modeling and methods that 
can take into account uncertainty in data and modeling, and 
incorporate biological and process constraints can provide a 
useful qualitative description of the responses of metabo- 
lism to metabolic engineering manipulations. The qualita- 
tive consistency of experimental results with trends sug- 
gested by optimization calculations using a {log)linear 
model indicate the utility of these mathematical tools for 
metabolic design. 

APPENDIX 

Details of the construction of the (log)linear model can be 
found elsewhere (Hatzimanikatis et al., 1996; Hatzimani- 
katis and Bailey, 1997a). Here, the structure and the param- 
eters of the (log)linear model used for the yeast glycolysis 
are summarized. The nonlinear model used to construct the 
(log)linear model can be found elsewhere (Hatzimanikatis 
and Bailey, 1997a; Schlosser et al, 1994). The mass bal- 
ances on the metabolites of the system may be written as: 



dz 

-—^Nez + NIIq 
at 



(1) 



and the metabolic outputs can be written as: 

w = 3ez + Hllq (2) 
where, z, q, and w, are the logarithmic deviations of the 




control 



PYKbs 



Figure 7. Specific glucose uptake (A) and cthnriol production rate (B) at 
2% initial glucose concentration of anaerobic, resting £. coli KO20 Strains 
with the plasmid p7Vr°°A only (control), over-expressing pjkF from E coli 
(PYKec), and overexpressing pyk from 8. sieamtkennopkilus (PYKbs). 
Different shading indicates different time points 'after the start of the ex- 
periment; 3 h (light grey), 5 h (grey), and 7 h ; (dark grey). Cells were 
harvested from the late exponential growth phase of aerobicaJly grown 
cultures. The indicated values represent the mean value of at least two 
independent experiments, and error bars indicate: the standard deviation. 

metabolite concentrations, the enzyme levels, and the meta- 
bolic outputs, respectively: 

Zi = ln(jt/X|,e) 

q k =\n(p k /p kt0 ) : (3) 
w, = ln( h t /h li0 ) 
and N> H, e, and II, are matrices, defined as: 



}. 



' »».|>o 



n = {ir M |ir yj , 



} 



(4) 



The metabolite vector, x, is defined for this model as: 
x r = [C ln , F6P, FdP, PEP, ATP) 
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where the italic symbols indicate the intracellular concen- 
trations of the corresponding metabolites. 
The reaction vector, v, is here: 

V 7 * = [V in , V HICt Vp OLt V PFKt V GROt V GAPO , V PYK. v A77wI 

The elasticity matrix e can be written as a sum of two 
matrices: 



e* + e' 



(5) 



where the elements in matrix e A correspond to the substrate 
elasticities of the enzymes, that is, the sensitivities of en- 
zyme activities with respect to concentrations of their sub- 
strates, and the elements of matrix e r correspond to the 
regulatory elasticities of the enzymes, that is, the sensitivi- 
ties of enzyme activities with respect to concentrations of 
regulatory metabolites. 

Changes of the elements in matrix e r from non-zero val- 
ues to zero* or vice-versa, correspond to modifications in the 
regulatory structure of the system. 

The stable steady state used as the reference state from 
which the linear model is constructed is; 

xj = (0.080200, 2.56666, 9.29469, 0-015318, 0.819456] 

and the values for the correspemding steady state rates are: 

vj = [27.366669 27.366669 4.374271 22.992396 
18.804933 4.187464 37.609866 20.486395] 

For the matrix N in Equation (1) we have; 

341.230 -341.230 0 0 
0 10.6624 -1.70426 -8.95809 
0 0 0 2.47371 

0 0 0 0 

0 0 0 0 

0 -33.3961 -5.33802 -28.0581 

0 0 0 0 

0 0 0 0 

-2.02319 -0.450522 0 0 
0 546.738 -2455.27 0 
O 10.2201 45-8961 -24.0 J 

The elasticity matrices at the reference steady state arc: 



N = 



0 0 
0.55174 0 
0 0.20134 
0 
0 
0 
0 



o 

0 
0 

0.41708 0 
0 0.72084 
0 0.054201 
0 0 
0 0 



0 
0 
0 
O 
0 
0 



0 
0 

0.05252 
0.03270 
O 

-0.06921 



0,35162 0.07157 
0 1.0 



— o -0.Q4689 


0 


0 0 


0 0 


0 


0 0 


0 0.93425 


0 


0 0 


0 0 


0 


0 -0.230352 


0 0 


0 


0 0 


0 0 


0 


0 0.16094 


0 0 


0.04261 0 0.02834 


_ 0 0 


0 


0 0 _ 



We consider the quantities of the enzymes catalyzing the 
first six reactions as the manipulated variables. Therefore, 
we can write for the manipulated variables the vector. 



-IV, 



fttjtr 



Vm,H*> V m,POlj Vm.rFK* 



m.CRO 



n = 



and 



and for the matrix II we have: 

10000006' 

0 1000000 

0 0 10 0 0 0 0 

0 0 0 10 0 0 0 

0 0 0 0 10 00 

0 0000 1 00 

0 0 000010 
_00000001 

Ethanol production has been considered to be propor- 
tional to the rate of the reaction catalyzed by pyruvate ki- 
nase. Therefore, the parameters in Equation (2) that corre- 
spond to Vpyf are: 



and e 4 and c r as given above. 

The synthesis problem was considered as in Hatzimani- 
katts et al. (1996). 
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